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We present experimental, numerical, and theoretical evidence for a new mode of antiferromagnetic 
dynamics in nanoparticles. Elastic neutron scattering experiments on 8 nm particles of hematite dis¬ 
play a loss of diffraction intensity with temperature, the intensity vanishing around 150 K. However, 
the signal from inelastic neutron scattering remains above that temperature, indicating a magnetic 
system in constant motion. In addition, the precession frequency of the inelastic magnetic signal 
shows an increase above 100 K. Numerical Langevin simulations of spin dynamics reproduce all 
measured neutron data and reveal that thermally activated spin canting gives rise to a new type of 
coherent magnetic precession mode. This ” rotor “ mode can be seen as a high-temperature version of 
superparamagnetism and is driven by exchange interactions between the two magnetic sublattices. 

The frequency of the rotor mode behaves in fair agreement with a simple analytical model, based 
on a high temperature approximation of the generally accepted Hamiltonian of the system. The 
extracted model parameters, as the magnetic interaction and the axial anisotropy, are in excellent 
agreement with results from Mossbauer spectroscopy. 


I. INTRODUCTION 

Magnetic nanoparticles have important applications in 
modern technology, including magnetic data storage, fer- 
rofluids, magnetic resonance imaging, biotechnology, and 
biomedicine^. In addition, the study of nanoparticle 
magnetism have elucidated many fundamental scientific 
phenomena, like thermally induced fluctuations, leading 
to magnetization reversal, known as superparamagnetism 
(SPM). The basic theory for SPM was derived by Neel 
and Brown&i, who found that the SPM relaxation time 
is given by an Arrhenius-like expression, 

T = nexp (j^r)’ (1) 

where r 0 is a typical attempt frequency and KV is the 
magnetic anisotropy barrier. Eq. © is supported by 
numerous experimental studies, see e.g. Refs. HI and y. 
Below temperatures where the SPM relaxation is impor¬ 
tant, the magnetic dynamics is dominated by a coherent 
uniform precession mode, which can be considered as a 
q = 0 spin wave. This mode gives rise to a linear decrease 
of the magnetization with increasing temperature^. 

The magnetic dynamics of antiferromagnetic (AFM) 
nanoparticles has attracted much attention, because it 
displays much richer intrinsic dynamics than that of fer¬ 
romagnetic nanoparticles^, and because complications 
due to dipolar interactions between particles can be ne¬ 
glected. In AFM nanoparticles with uniaxial anisotropy, 
the uniform precession mode resembles a q = 0 AFM 
spin wave, and its frequency is influenced both by the 
anisotropy and by the strength of the magnetic exchange 


interaction^. 

The technique of inelastic neutron scattering have ear¬ 
lier been found to be efficient to investigate the mag¬ 
netic dynamics of nanoparticles, see e.g. Refs. US In 
particular, in 15 nm hematite (a-Fe 203 ) particles, the 
frequency of the uniform modes was found to decrease 
with increasing temperature. This was ascribed to an- 
harmonicities in the uniaxial anisotropy potential^. 

In this work, we present neutron diffraction data on 
8 nm hematite particles, which show an unexpected van¬ 
ishing of static magnetic order at temperatures as low 
as 150 K, despite the fact that the Neel temperature of 
bulk hematite exceeds 900 K. We analyse this result in 
the light of inelastic neutron data, which show a clear sig¬ 
nal from the uniform magnetic mode at all temperatures 
studied (up to 300 K) and an increase with temperature 
of the uniform mode frequency. We further investigate 
the physics of this system by a theoretical model and 
by numerical simulations, and we find that they lead to¬ 
wards the existence of a novel dynamic ” rotor “ mode in 
AFM nanoparticles. 


II. THE NANOPARTICLE SAMPLE 

Hematite is a common AFM mineral occuring fre¬ 
quently in soils as nanoparticles. The magnetic ions in 
hematite are Fe 3+ (s = 5/2), which are essentially ar¬ 
ranged in a hexagonal structure. The spins are ferro- 
magnetically ordered in the (a, b) planes, with alternating 
directions along the c direction. In bulk hematite above 
the Morin transition, Tm = 263 K, the spins are confined 
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to the plane by a strong planar anisotropy, while below 
7 m, the spins turn to point along the c-axis— 

In particles smaller than ~ 20 nm, the Morin transi¬ 
tion is suppressed. At low temperatures, the spins align 
along an easy direction within the (a, b) plane due to 
a weak axial anisotropy, presumably from the particle 
surface. The sublattices display a small canting from an¬ 
tiparallel within the plane (in bulk ^0.1 o1 —), due to the 
Dzyaloshinsky-Moria interaction. Our simulations show 
that this small term has little effect on the properties 
we investigate, and it is generally irrelevant in the fol¬ 
lowing discussion. A later publication will discuss this 
approximation in more detail^. 

The nanoparticle sample of hematite was characterized 
earlier—. The particles are phosphate coated at the sur¬ 
face. Electron microscopy and X-ray diffraction show 
that the particles have an average diameter of 8 nm. 
The phosphate layer is very thin, most likely a mono- 
layer, but sufficient to minimize magnetic interparticle 
interactions^. 


III. NEUTRON SCATTERING EXPERIMENTS 

Neutron diffraction experiment were performed at the 
Paul Scherrer Institute (CH), using the two-axis pow¬ 
der diffractometer DMC and the triple-axis spectrometer 
RITA-2. The DMC experiments were performed with a 
wavelength of 4.2 A (4.7 meV), giving an effective energy 
resolution wider than 5 meV (FWHM), while the RITA-2 
experiments were performed at 4.7 A (3.7 meV) with an 
energy resolution of 0.12 meV (FWHM). 

The inelastic neutron scattering experiments were per¬ 
formed with RITA-2 in the monochromatic imaging 
mod o 21 ’ 22 . These experiments were performed with a 
constant final energy of 3.7 meV, also here with energy 
resolution of 0.12 meV (FWHM). 

Quasielastic neutron experiments were performed at 
the backscattering spectrometer BSS at Research Center 
Jiilich, using a neutron wavelength of 6.271 A and an en¬ 
ergy resolution of 1.0 /reV (FWHM). Of the 14 detectors 
of BSS, we used number 2 to 5, corresponding to q -values 
of 1.9, 1.7, 1.5, and 1.3 A -1 , respectively. Each detector 
covers a g-range of around 0.2 A -1 . Detectors 3, 4, and 5 
correspond to the structural (102) peak and the magnetic 
(101) and (003) peaks, respectively. 


IV. NEUTRON SCATTERING DATA 

Fig. [I] shows the results of the diffraction experiments. 
The AFM hematite peaks (003) and (101) are clearly vis¬ 
ible at q -values of 1.37 A -1 and 1.51 A” 1 , respectively, in 
both data sets. However, the temperature dependence of 
the magnetic scattering intensity is most different in the 
two experiments. In the two-axis experiment the peak in¬ 
tensity at 300 K is reduced by 10% relative to the value 
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FIG. 1. (color online) Neutron diffraction data around the 
magnetic (003) and (101) AFM Bragg peaks, corresponding 
to q = 1.37 A -1 and q = 1.51 A -1 , respectively, taken at 10 K 
and 300 K. (top) data from the two-axis powder diffractome¬ 
ter DMC. (bottom) data from the triple-axis spectrometer 
RITA-2, using energy analysis of the scattered neutrons. 

at 10 K, whereas the triple-axis experiment shows a re¬ 
duction of 90%. The only difference between these two 
experiments is the energy resolution, indicating that the 
hematite spins become almost completely dynamic with 
an energy scale in the range between the respective en¬ 
ergy resolutions of the two instruments. 

These initial observations prompted us to re-analyse 
previously published inelastic neutron scattering data 
from the same sample in the temperature range 10- 
200 K—, as well as unpublished data up to 300 K. Fig. [2] 
shows our data, displaying the signal at a scattering vec¬ 
tor of q = 1.37 A -1 , corresponding to the AFM (003) 
reflection. The strongest feature is the central elastic 
peak that stems from elastic incoherent and quasielastic 
background and a narrow quasielastic SPM signal. 

The uniform modes are seen in the data as broad side 
peaks at \huj a \ ~ 0.3 meV. We model these peaks by 
a damped harmonic oscillator (DHO) functions!^. Be¬ 
fore correcting for experimental resolution and constant 
background, the model reads 

J(w) = Aj nc <5(u;) + C (2) 

Aq e rqe/7r A SPM r/7r 2A DH o-P(a;)ia 2 7 a /7r 

r 2 e + W 2 r 2 + w 2 (w 2 - w 2 ) 2 + 47 2 w 2 ’ 

where Tiui is the neutron energy transfer, A; nc , A qe , 
Aspm, and A^ho are the integrated intensities of the in¬ 
coherent elastic, the quesielastic incoherent (from mobile 







3 



hui [meV] 



hui [meV] 


FIG. 2. (color online) Inelastic neutron scattering data at 
q — 1.37 A' 1 , corresponding to the (003) AFM Bragg peak, 
measured at 50 K (top) and 150 K (bottom). Counts are 
normalized to monitor and plotted on logarithmic axis vs. 
energy transfer, hco. The solid red lines show the best fit 
to the model m convoluted with experimental broadening. 
Individual components of the fits are shown by dot-dashed 
black lines, solid green lines, and dashed blue lines, which 
represent the incoherent elastic background, the quasielastic 
SPM signal, and the DHO component, respectively. 


H 2 0), the quasielastic magnetic (from SPM), and the 
magnetic DHO components, respectively. r qe , T, and 
are the widths (HWHM) of the quasielastic incoherent, 
the quasielastic magnetic, and the inelastic peaks, respec¬ 
tively. D(w) = w[n(w) + l\/[k&T\ is the detailed balance 
factor—, with n(w) being the Bose factor and fce being 
the Boltzmann constant, as detailed ini^. The values of 
Ainc, A qe , and r qe are found from constant-g scans at 
non-magnetic q- values, leaving seven free fitting param¬ 
eters for the magnetic scattering: AgpM, T, Adhcx 
7 a , C and a possible offset of the elastic line due to im¬ 
precision in alignment. 

Experimental data taken at 10 K have intrinsically low 
signal-to-noise ratio. The parameter values for the in¬ 
elastic signal then becomes unreliable, since the fit tends 
to model non-Gaussian tails in the resolution function 
for the large elastic signal, rather than the tiny inelastic 
magnetic signal. It was therefore necessary to disregard 
data taken at this temperature in the data analysis. 

The quality of the fit of the model (0 to data is cru¬ 
cial to this work. Some correlation could be expected 
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FIG. 3. (color online) The quality of the model fit to the 
inelastic data, measured as the unreduced chi-square value, 
Xunred i shown for the parameter plane spanned by and 
7 a. At each point in this plane, the value of the five other 
free fitting parameters have been fitted to the 115 data points. 
The 95% confidence interval, corresponding to AxSnred = 2-3, 
is marked by the solid black line. 




FIG. 4. (color online) Position (u a ) (top), and width ( 7 ^) 
(bottom) of the inelastic DHO signal as a function of temper¬ 
ature. Black circles are measured values and blue diamonds 
are results from the Langevin simulations, described in the 
text. Error bars on the simulated data points are comparable 
to the symbol sizes. In the top panel, the dot-dashed green 
curve in the hui a plot is the previously published result for 
the uniform mode^, the dashed red line is the more accurate 
version of the same calculation—, and the solid black curve is 
our analytical result for the ”rotor" mode, eq. GD- 
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FIG. 5. (Color online) quasielastic neutron scattering data 
from BSS, Jiilich, taken close to the magnetic (003) AFM 
Bragg peak, (top) raw data taken at T = 50 K with a 
Lorentzian model fit to data as explained in the text, (bot¬ 
tom) Fitted Lorentzian line widths of the superparamagnetic 
signal as a function of temperature. The solid lines are fits to 
the Neel-Brown law as described in the text. 


between the free parameters. For the important DHO 
signal, two parameters are of the same energy scale of 
a few tenths of a meV: 7 a and ui a , and we have there¬ 
fore investigated how these two parameters correlate by 
fixing their values and performing a fit of the other five 
parameters. Then, we record the goodness of fit given 
by the value of the unreduced chi-square function, X„ nred 
function. We have done this for a range of parameters, 
spanning the (j a , w a )-plane with the results presented in 
Fig. [3] The confidence interval in such a figure is given 
by the condition^ ^X'unred < 2.3 and is marked by the 
black solid line in the figure. We see that the values of 
the two parameters are essentially uncorrelated and that 
the fitted parameters are robust; except that the 7 a pa¬ 
rameter may have a somewhat asymmetric errorbar. 

Our central findings of experiments and model fit, the 
energy of the uniform mode hui a , is shown in fig. [I] oj a 
is found to increase with temperature, in direct contrast 
to the results from 15 nrn particles^. 

The width of the central peak was further studied 
by high-resolution quesielastic measurements around the 
magnetic (003) peak, addressing superparamagnetism 
fluctuations. One example of the raw data is presented 
in fig. [5] (top) together with a model fit that contains 
a Lorentzian line shape convoluted with the resolution 


function, plus a constant background. In fig. 0 (bottom), 
the Lorentzian linewidths are plotted as a function of 
temperature. The width is seen to increase with tem¬ 
perature up to around 150 K, where the magnetic signal 
becomes wider than the instrument energy window and 
no useful information can be extracted. The data up to 
150 K is well modeled by the Neel-Brown law ©, giv¬ 
ing an anisotropy value of KV/ks = 194(40) K with 
r 0 = 9(2) x 10 " 12 s. 

The remainder of this article is devoted to the expla¬ 
nation and modeling of the two unforeseen results: The 
vanishing of the diffraction peak at room temperature in 
the triple-axis experiments and the increase of precession 
frequency with temperature. In the discussions, we will 
also include the quasielastic neutron scattering data. 


V. ANALYTICAL THEORY 

A derivation for the spectrum of a fluctuating ferro¬ 
magnetic nanoparticle in a uniaxial anisotropy was ear¬ 
lier derived by Wiirger—. However, in this scenario, the 
frequency of the uniform modes decrease with increasing 
temperature, in contrary to our observations. For this 
reason, we decided to develop an analytical understand¬ 
ing for this AFM system with two anisotropies. 

The generally accepted microscopic magnetic Hamilto¬ 
nian for hematite nanoparticles reads^: 

n = E • s ; - *1 E( s i ) 2 - ^ E«) 2 > ( 3 ) 

ij i i 

where i and j are atomic indices, k± < 0 is a large pla¬ 
nar anisotropy, and K 2 > 0 is a smaller axial anisotropy, 
which is linked to the anisotropy barrier in Eq. © 
through KV/k-Q = K 2 Ss'N, where s' = s — 1/2— J.y 
is the exchange coupling constants. In this notation, y 
and ^ lie within the basal hematite (a, b) plane, while x 
is along the c-axis. 

To describe the low-energy dynamics, we follow 
Wiirger— and assume collective motion (q = 0) of the 
spins in each of the two sublattices (A and B). This is 
justified by the fact that the highest temperature in our 
study, 300 K, is far below the Neel temperature of the 
system and therefore the spin waves have reduced the 
ordered moment of the sublattices only slightly. In addi¬ 
tion, the finite size of the system results in a quantization 
of the spin wave spectrum, and due to the steepness of 
the hematite spin wave dispersion, only few spin waves 
states are in fact allowed below energies corresponding 
to 300 K21. 

We define the number of spins in the two sublattices 
as, Na = Nb and define the uncompensated fraction £ = 
Na/Nb■ We initially consider the case £ = 1. The spin 
dynamics can be fully described by the two superspins 

£> A = ^E Si , S B = 

ieA ieB 


( 4 ) 
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In terms of Sa and Sb, the Hamiltonian becomes 
H = J'S a -Sb-k' 1 (S¥ + S% 2 ) 

-k' 2 (S z a 2 + S Z B 2 ) - 7 B-(Sa + Sb), (5) 

To simplify, we define Sa = Sa /S, and similarly for Sb- 
We then reach 

n = 7 ^[-B ex S A • Sb + + S£ 2 ) 

-^{Sx + ^ 2 )l ( 6 ) 

where the gyromagnetic ratio is 7 = gg, B /h and H ex 
B\ > B 2 . Here, Hi and B 2 are the planar and uniaxial 
anisotropy fields, respectively, defined by gg B Bi = 2\ni\s' 
for i = 1,2. The exchange field is gg B B ex = 2 zJ nn s, 
where the factor of 2 is due to double-counting. 

In bulk hematite, H ex = 827 T 12 i 19 i 26 , From the neu¬ 
tron scattering data presented in Fig. [2] and the model 
Eq. [TO] we find B -2 = 4.8 mT. B\ cannot be determined 
from the present neutron scattering data and we there¬ 
fore use the value Hi (8 nm) = Hi (16 nm)i£, Hi « 46 mT 
at low temperatures, gradually increasing to ~ 66 mT at 
300 K. 

The exchange term is minimal when Sa and Sb point 
in opposite directions; the Hi anisotropy is minimized 
if they lie in the (y, z) plane, and the B -2 anisotropy will 
attain its minimum when Sa and Sb are parallel to the pr¬ 
axis. Hence we have two minimum energy configurations 
S A = (0, 0, ±1), S B = (0, 0, =Fl). 

Since the magnitudes of the superspins are much larger 
than h, quantum mechanics plays no role, and the dy¬ 
namics is governed by the classical equations of motion: 

d 9® 

-^ = [S\U], (7) 

where the right hand side is a Poisson bracket, obeying 
the rule [ST S v ] = S z , etc. Using H from f 6 ]), we get 

^ = 7 SaxBa,5=7SbxB b , ( 8 ) 

dt at 

where Ba is the resulting field on sublattice A originating 
partly from the anisotropies, partly from the magnetisa¬ 
tion of sublattice B: 

B A = - H ex S B + B 1 S A e x + B 2 S z A e z , 

B b = - H ex S A + HiHge* + B 2 S^e z , (9) 


These are the two uniform modes calculated ear¬ 
lier, and also observed directly with inelastic neutron 
scattering* ^! .! ! .. The model can be extended to include 
£ ^ 1, see Ref. 0. 

The temperature dependence of the excitation energies 
at low temperatures is given by 

u a ,p(T)=u a ,p(O)0 z ), (11) 

where the ordered sublattice moment ( S z ) decreases 
monotonically with temperature, according to Boltz¬ 
mann statistic :) 12 ' 17 . 

It turns out that the system supports another approx¬ 
imate mode, which we will denote the “rotor” mode, de¬ 
scribed in the following. We assume that the magnetic 
moments move slightly out of the (y, z) plane, but not 
with strictly antiparallel sublattices, e.g. 

S A = (9 A , S v , S z ) S B = (0b, -S v , -S z ). (12) 

We first note that the canting angles described here are 
out of the basal plane, and therefore unrelated to the in¬ 
plane canting angle caused by the DM-interaction. With 
our ansatz , the equations of motion become 

d0A R Q Q 

— = 7 B 2 S y S z 

dfe nee 

— = 7 B 2 S y S z 

1 Q 

-TT = -lB ex {0 A + 0 B )S Z 
dt 

^- = jB ex (e A + e B )s y , ( 13 ) 

where we neglect terms of order Hi in comparison to 
terms of order H ex . Assuming that the angles d A and d B 
are constant in time, the magnetic moments will perform 
full 27 t rotations around the cc-axis with a frequency given 
by 

Wrot = lB ex (6 A +0 B ). (14) 

This rotation is illustrated in Fig. [ 6 ] for the case 0 A = 0 B . 
We denote the sublattice canting angle 8 = \9 A +9 B \- The 
9 dependence of the overall system energy stems from the 
exchange interaction, Eg oc JS 2 cos (8). From Boltzmann 
statistics, we can now calculate the thermal average value 
of 9 2 to 

(8 2 ) « 2fc B T/( 7 HH ex ). (15) 


The scalar S\ is the projection of S A onto the x-axis. 
Considering only motion with small deviations from the 
minimum energy configuration, we obtain two normal 
modes with frequencies^ 


uj a — 7 \J (2H ex + Hi + B 2 )B 2 ~ 7v / 2H ex H2 (10) 
W/3 = 7 \J (2Hex + H 2 )(Hi + B 2 ) ~ 7\/2Hex(Hi + B 2 ). 


We next investigate the assumption that the relative 
rate of change of the angles, 8/9 , is small compared to 
w ro t- Neglecting thermal fluctuations we obtain 


M~_^~ 5x10 -3^ 

Wrot 2H ex 0 2 T[K] ’ 


( 16 ) 


where we have used the values of B 2 and H ex given in 
the following. For nanoparticles corresponding to 8 nm 
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FIG. 6. (Color online) Illustration of two magnetic modes in 
an antiferromagnetic nanoparticle with a vertical easy axis. 
The A (B) sublattice points in the +z (—z ) direction and de¬ 
scribes a blue (red) trajectory: a) a uniform in-plan magnetic 
mode; the ooamode; b) the new ’’rotor mode". 


hematite, N A ~ 5000, and for T ~ 100 K, (OH gives 
0/{0u TOt ) = 0.3 < 1. Hence, 9a and 9 B maintain their 
values for sufficiently long times that the rotor mode is 
significant. 


VI. NUMERICAL SIMULATIONS AND 
RESULTS 


To veryfy the soundness of the analytical approxi¬ 
mations, and to obtain additional detail in the model, 
we numerically simulate the sublattice dynamics, using 
Langevin dynamics. Here, the equation of motion is a 
stochastic Landau-Lifschitz equation, very similar to the 
method applied earlier to describe superparamagnetism 
in ferromagnetic nanoparticles— In the numerical cal¬ 
culations, we have lifted the requirement £ = 1 and also 
included the Dzyaloshinsky-Moriya interaction. In this 
scheme, the equations of motion for the two sublattices 


become: 

^-S A = 7 Sa x (B a + b) - A| 7 |S a x (Sa x B a )(17) 
^S B = 7 S b x (B b + b) - A| 7 |S b x (S b x B b ), 
where 

B a = £B ex S B + rhS x A e x + 

B 2 S z A e z - £B D (£|e" - S y B e z ), 

B b = Bex S A + + 

B 2 S z B e z + B D (S z A e y - S y A e z ) (18) 

Bo is the Dzyaloshinsky-Moriya field, where we use the 
bulk hematite value of Bp = 2.1 T— . In (1171) . b is a small 
random field representing thermal fluctuations, which are 
isotropic and uncorrelated in time. The term A repre¬ 
sents dissipation. The fluctuation-dissipation theorem-22 
is used to determine in absolute units the simulated tem¬ 
perature: 


* b T = | 7 | guBS ^2^ • ( 19 ) 

Note that the particle size is present here through N A . 
The fluctuating field is uncorrelated in time, {b(t)b(t')) oc 

The equations of motion m were integrated numer¬ 
ically using the Runge-Kutta method of second order— 
with a time step of 1 fs (1000 THz), much faster than 
typical frequencies of the uniform mode io a (0.5 THz). 
The method was tested against analytical results in the 
limit of vanishing fluctuation and dissipation^^, and for 
the one-sublattice (ferromagnetic) cas e 2 - 2 i— . 

Examples of raw simulation data for the two-sublattice 
case are shown in Fig. [7] The motions of the sublat¬ 
tice magnetizations were, in turn, Fourier transformed 
to obtain the power spectra, which through the Wiener- 
Khintzine formula— yield the simulated scattering inten¬ 
sities. Since the scattering vector lies in the ^-direction, 
we sum the power spectra of the y and z spin compo¬ 
nents to model experimental data. The simulated data 
was convoluted with the experimental resolution and fit¬ 
ted in the same way as the experimental data, using m, 
resulting in good agreements over three orders of mag¬ 
nitude; see Fig. [ 8 ] The deviations found at the small- 
amplitude data at | hut\ > 0.8 meV are irrelevant for our 
conclusions and we neglect them in the following. 

In the simulations we have adjusted four parameters: 
The uniaxial anisotropy field B 2 , the exchange field B ex , 
the uncompensated moment fraction, £ and the dissipa¬ 
tion parameter A to find the values that yield the best 
agreement between simulations and experiment with re¬ 
spect to the temperatures dependences of 7a and Wq— . 
We find B 2 = 9.5 mT, corresponding to KV/kn = 169 K. 
Bx and £ cannot be determined independently: any 
value of Bx between 600 and 1000 T has a correspond¬ 
ing value of £ where simulations and data match. We 
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FIG. 7. (color online) Raw data from the Langevin simula¬ 
tions, showing the time development of the three spin com¬ 
ponents at temperatures of 40 K (top), 150 K (center) and 
300 K (bottom). 9(t) represents the canting and is displayed 
in degrees. The encircled areas in the center and bottom panel 
show fast coherent oscillations in the j/s-plane correlated with 
a large value of the canting angle. 


therefore fix the exchange field to the bulk value, Bx = 
827 T, and thus find £ = 1.010. Finally, we determine 
A = 6 x 10~ 4 . 


VII. DISCUSSION 
A. Existence of the rotor mode 

We first discuss our hypothesis of the existence of a 
rotor mode by comparing data from neutron scattering 
experiments, numerical simulations, and analytical cal¬ 
culations. 

Fig. [4] (top) displays the temperature dependence of 
the mean precession frequency, w Q , defined through ©• 
The agreement on the temperature behaviour of w a be¬ 
tween experiment and simulation is excellent. 

The value of the DHO peak width 7 a , displayed in 
Fig. 2] (bottom), is generally higher in the experimental 
data than in the simulations. We ascribe the higher ex¬ 
perimental value of j a to a distribution of particle sizes 
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FIG. 8. (color online) Power spectra from the Langevin sim¬ 
ulations at 150 K. Top: The individual components, S x (solid 
red), S y (dashed blue), and S z (dash-dotted green). All three 
oscillator modes are clearly visible, as indicated in the figure. 
Bottom: Sum of the y and a components (solid black) with a 
fit (solid red). The individual components of the fit are the 
SPM and the DHO signals, given by the dash-dotted green 
and dashed black lines, respectively. 


and shapes, which results in a spread in uj a . This would 
appear in the fit as an overall increase of ") a , most promi¬ 
nent at the lowest values of this parameter, as the parti¬ 
cle size distribution effect is presumably constant and the 
width of the two broadenings should add in quadrature. 
This is in good overall agreement with the experimental 
observations. 

The simulated values of u a are in excellent agree¬ 
ment with the neutron scattering data. At low tempera¬ 
tures, T < 50 K, the analytical uniform mode prediction 
agrees well with the simulations. At higher temperatures, 
the uniform mode prediction strongly underestimates the 
value of u) a , while our analytical results for the rotor 
mode overestimates this value. We believe the reason for 
this to be the following: Both the uniform mode and the 
rotor mode give rise to a DHO signal, which due to the 
relatively short life time of the excitations are both quite 
broad. Experimentally we cannot resolve the two indi¬ 
vidual DHOs, and a single DHO has therefore been used 
to fit the data. The observed frequency is thus a weighted 
average of the frequency of the rotor mode, w ro t and that 
of the uniform mode, ui a , in general agreement with the 
experimental and simulated data. 

Additional insight is found by inspecting the simulated 





































































































































time evolution of the magnetic dynamics, shown in Fig. [7] 
The 40 K data show low-amplitude, high-frequency os¬ 
cillations in S x and high-amplitude, low-frequency os¬ 
cillations in S v , representing the two uniform modes 
m- The easy-axis component, S z , displays nearly con¬ 
stant values, interrupted by sudden (SPM) magnetiza¬ 
tion reversals. These three observations are in accor¬ 
dance with the present understanding of magnetic dy¬ 
namics of nanoparticles. The 150 K data contains peri¬ 
ods with much less regular motion of the sublattice spins. 
In particular, both S z and S y show oscillations. The 
highlighted regions show time intervals with coupled fast 
oscillations in these parameters, correlated with a high 
value of 9. These fast oscillations can be seen as re¬ 
peated, coherent SPM relaxations and rarely appear at 
the lowest temperatures. At 150 K they appear at least 
once per simulated nanosecond, and dominate at larger 
temperatures, e.g. at 300 K. These events are indeed 
occurences of short-lived rotor dynamics. 

Our explanations are supported by the power spectra 
shown in Fig. [HJ where the SPM behaviour is found in 
spin components along the easy 2 -direction, the DHO 
signals are seen in the y and x components, while the 
rotor mode is seen as a high-frequency shoulder in the z 
signal. The expected rotor contribution to the spectrum 
of S x is not clearly visible due to the broadening of the 
signal from the uj a mode. 


B. Comparison with earlier work 

In earlier work by our group, we found a value for the 
energy barrier KV/kB = 250(30) K using the ratio be¬ 
tween elastic and inelastic signal in the neutron scattering 
data from the same sample, and KV/k b = 335(35) K us¬ 
ing Mossbauer spectroscopy data of the same sample^. 
However, this work used an effective 3D model, in which 
the in-plane anisotropy is neglected, leading to an overes¬ 
timation of KV. Assuming that the in-plane anisotropy 
is identical to the value for 16 nrn particles^, and that 
£ = l,the values become KV/kB = 148(18) K from 
the neutron data and I\V/kB = 211(20) K from the 
Mossbauer data. From the values of u> a , we furthermore 
found an uncompensated moment of £ = 1.010 using the 
uniaxial model. 

From the fit of the backscattering data to the Neel- 
Brown model, we find in the present work a value of 
the energy barrier, KV of 194(40) K, in agreement with 
the Mossbauer data. In addition we obtain the prefactor 
value Tq™ = 9(2) x 10 -12 s. 

Tuning the simulations to obtain best agreement be¬ 
tween the measured and simulated values of hco a and , 
we reach a ” modeled refinement” of the neutron data, gv- 
ing KV/kB = 169 K and £ = 1.010, in very good agree¬ 
ment with earlier work. The simulated backscattering 
data using the same simulation parameters also follow the 
Neel-Brown law very closely, leading to KV/kB = 169 K, 
but with a slightly higher prefactor value, Tq™ = 6.8 x 


10~ 12 s. Hence, our data and simulations show a very 
good internal consistency, except for the r 0 value. This 
relatively minor difference, we ascribe to anharmonicities 
in the anisotropy potential, not included in the simula¬ 
tions. 


C. Uncompensated moment and coupling of modes 

Contrary to expectations from simple theory, our sim¬ 
ulations show coupling between the Lo a and up modes. 
This is most obvious in the power spectra in Fig. [51 where 
the x-component has a broad central peak. At lower tem¬ 
peratures, this signal consists of two peaks at ±hui a ; the 
peaks gradually merge as the temperature is increased. 
By varying the simulation parameters, we found the cou¬ 
pling to be due to the uncompensated moment. At £ = 1, 
the modes seem to be completely uncoupled as expected. 
However, our simulations with £ = 1.010 in our view rep¬ 
resents closer the physical reality for a typical nanoparti¬ 
cle, and it is satisfying that even with this complication 
we are able to reproduce the experimental observations 
so accurately. 

D. Absence of elastic scattering 

As a final test of the rotor model, we return to the ini¬ 
tial observation of the disappearance of the elastic neu¬ 
tron scattering signal when measuring with energy anal¬ 
ysis on a triple-axis instrument. 

The top panel of Fig. [9] shows the elastic magnetic sig¬ 
nal when integrated over the energy resolution of BSS, 
RITA-II and DMC, respectively. With the resolution of 
BSS of ~ 1 /xeV, there is almost no elastic signal at tem¬ 
peratures above ~ 100 K. At RITA, the elastic signal at 
300 K has decreased to approximately 20% of the low 
temperature signal, while within the DMC resolution, 
most of the signal is still seen as elastic at 300 K. 

Overlayed on this, Fig. [9] shows the corresponding in¬ 
tensities from the numerical simulations, scaled only by a 
common factor. The agreement between experiment and 
simulation is striking for the DMC and RITA-2 data, 
even though A§pm was not used in the tuning of the 
simulations. This gives additional confidence in the sim¬ 
ulated model. In the backscattering data, the experimen¬ 
tal value of the elastic signal is in general higher than the 
simulations. This may well be a combination of two ef¬ 
fects: a) the difference in r 0 value, b) a distribution of 
particle sizes, and therefore anisotropy barriers. Neither 
of these effects are included in the numerical simulations. 

The bottom panel of Fig. [9] shows the total area of the 
fitted quasielastic magnetic signal, Aspm as a function of 
temperature. This value has been corrected for inelastic 
scattering and should thus be equal for both the BSS 
and RITA experiment, which is indeed the case within 
the uncertainties. Agreement with the simulations is also 
found in this case. 
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FIG. 9. (color online) Temperature dependence of the elastic 
part of the inelastic neutron scattering signal, (top) Back¬ 
ground subtracted diffraction data from the backscattering 
spectrometer BSS, the triple-axis spectrometer RITA-2 and 
the diffractometer DMC. (bottom) The area of the quasielas¬ 
tic superparamagnetic signal as obtained from modeling the 
data from BSS and RITA-2. 


servation. 

We speculate that the rotor mode may be present in 
other nanoparticles with a strong planar anisotropy, e.g. 
NiO, although the temperature effect in this system was 
found to be less pronounced^. 


VIII. CONCLUSION 

We observe an increase with temperature of the mean 
frequency of the uniform mode in weakly interacting an¬ 
tiferromagnetic 8 mu hematite particles. Simultaneously, 
we observe a disappearance of the strictly elastic signal 
from static magnetic order. Numerical Langevin simula¬ 
tions of a fairly detailed model is able to reproduce the 
features of the experimental data with sufficient accu¬ 
racy. In combination with an analytical model, the sim¬ 
ulation show that both observations can be understood 
by the presence of a new fast dynamical mode - the “ro¬ 
tor” mode, which can be seen as a coherent superpara¬ 
magnetic relaxation. In contrast to other uniform modes 
in magnetic nanoparticles, the rotor mode is driven by 
temperature-induced canting of the sublattice spins and 
causes a periodic magnetisation reversal that transforms 
the spin system to become purely dynamical. Our results 
are in good quantitative agreement with earlier neutron 
scattering and Mossbauer data. 

We predict this mode to be of general relevance for 
antiferromagnetic nanoparticles with strong easy-plane 
anisotropy. 


E. Perspectives of the rotor mode 

Our analytical and numerical results provide a qualita¬ 
tive explanation of the experimentally observed increase 
of u> a : When thermal excitations cause a canting of the 
sublattice spins out of the basal plane, a torque from the 
exchange interaction will appear. This dominates the 
axial anisotropy, B 2 , and causes the spins to perform co¬ 
herent precessions within the easy plane, hence inducing 
periodic reversals of S z . The frequency of the preces¬ 
sion depends on the canting, and hence on temperature. 
This displaces the neutron scattering signal from (quasi-) 
elastic to inelastic, in accordance with experimental ob- 
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